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Recent Results on Coding and Nonlinear 
Equalization at the University of Arizona 


V 

w 

w 


William Ryan 
University of Arizona 


March 8, 2001 
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Outline 


• early turbo code designs (serial and parallel) 

• more recent turbo code designs: double turbo DPSK 

• initial work on low-density parity-check codes 

• combined nonlinear equalization and decoding 



Bandwidth-Efficient Parallel and Serial 
Concatenated Codes for QPSK Modulation: 

Consideration for the CCSDS Near-Earth Standard 

William E. Ryan, University of Arizona 
Omer Acikel, New Mexico State University 

March 4, 2001 


OUTLINE 


• Bandwidth-efficient signaling with MPSK 

• The proposed codes - parallel and serial concatenations 

• Performance results for the codes 
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BANDWIDTH-EFFICIENT SIGNALING 
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• bandwidth-efficiency for coded QPSK (r = code rate): 

£qpsk, coded 2r (bits/sec) /Hz 
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• bandwidth-efficiency for coded 8PSK (r = code rate): 

Zspsk, coded = 3 r (bits/sec) /Hz 


Example 

• rate 2/3 trellis coded 8PSK: 

£s psk, coded = 2 (bits/sec) /Hz 

• rate 15/16 turbo coded QPSK: 

£qpsk, coded = 1-875 (bits/sec)/Hz 

• but, the coded QPSK scheme is more power efficient and has 
a more robust receiver 
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THE PROPOSED CODES 
The Parallel Concatenated Convolutional Codes 


AWGN 



Vi 




Diagram of PCCC encoder. 


• we consider rate 3/4, 7/8, and 15/16 PCCCs, with memory m = 
3 component recursive systematic convolutional codes 
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The Serial Concatenated Convolutional Codes 
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SCCC encoder with differential encoder as rate 1 inner code. 



SCCC encoder with two rate 1/2 RSC encoders. 
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(i) Top Diagram: a punctured outer rate r 0 = k/(k + 1) RSC 
encoder and a 7* = 1 (2-state) differential encoder 


( 2 ) Bottom Diagram: two m = 3 RSC encoders, with the inner 
and outer RSCs punctured to achieve rates r{ = 2k/(2k+l) and 
r 0 = (2 k + l)/(2 k + 2), respectively 

• We consider rate 3/4, 7/8, and 15/16 SCCCs. 
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PERFORMANCE RESULTS 
PCCC Performance 

• decoder: iterative BCJR-APP algorithm with 15 iterations 

• information word size: 10,000 bits 

• code rates: 3/4, 7/8, and 15/16 
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Rate 3/4 PCCC bit error rate performance. 
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Pb (BER) 


w 

V 

w 

w 


w 

W 

W 




V 


V 



Rate 7/8 PCCC bit error rate performance. 
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Pb (BER) 
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SCCC Performance 

• decoder: iterative BCJR-SISO algorithm with 15 iterations 

• information word size: 10,000 bits 

• code rates: 3/4, 7/8, and 15/16 
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Rate 3/4 SCCC bit error rate performance. 
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Eb/No (dB) 


Rate 7/8 SCCC bit error rate performance. 
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Pb (BER) 



Rate 15/16 SCCC bit error rate performance. 
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Performance Comparisons: PCCC, SCCC, Convolu- 
tional, and Trellis-Coded 8-PSK 



Comparision of rate 3/4 codes: PCCC, 8/2 SCCC, and a memor 3 r 6 

convolutional code. 
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Pb <BER) 



Comparison of rate 7/8 codes: PCCC, 8/2 SCCC, and a memory 6 

convolutional code. 
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8/2 SCCC 
Convolutional, m=6 
r=2/3 TC 8-PSK 
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Comparision of the rate 15/16 PCCC and 8/2 SCCC codes, and 

rate 2/3 trellis-coded 8PSK. 
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C and r (bits/channel use) 


Comparison to QPSK Capacity Curve 


• the Eb/No required to achieve a bit error rate of 10 5 is plotted 

« r="'i £ 

against the QPSK capacity curve 



Bandwidth-/Power-Efiiciency operating points for the codes 


presented in this paper. 
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Convolutional Double Accumulate Codes 
(or Double Turbo DPSK) 

Rajeev Ramamurthy and William E. Ryan 
ECE Department, University of Arizona 

November 15, 2000 




n = Nt 0 hit pseudo-randbm irtedeaver II " = Nr 0 bit pseudo-random de-interleaver 

- Outer convoiudonal code (rate r 0 < 1) D 0 = CXrter APP decoder. 

Ej! , Inner axuirulate codes (rale r t l ,% = !) D tl , D l2 = Inner AFP decoders. 


Encoder and decoder for double turbo DPSK. 

























Eb/No (dB) 

FER’s of rate serial, parallel codes 





Initial Results on 
Low-Density Parity-Check Codes 


Michael Yang and William Ryan 
University of Arizona 

March 2001 



low-density parity-check codes were invented by Gallager 
circa 1960 at MIT for his doctoral dissertation 

approximately five years ago, they were independently re- 
invented by MacKay (Cambridge University) and others 

the idea behind these codes (we assume binary codes): 

1. recall all codewords in a linear code form the 
nullspace of a so-called parity-check matrix H: 

Hc t = 0 whenever c is a codeword 

2. we say the binary matrix H is low-density if it is 
sparse, that is, if it contains only a small fraction of 
ones 

3. the low-density property leads to very powerful 
codes with decoders lower in complexity than turbo 
decoders 
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• the initial designs on LDPC codes (Gallager and MacKay) 
involved random designs with certain row and column 
constraints (e.g, # l’s per column k c > 2 and the # l’s per row 

K > k c ) 

• the problem with randomly generated codes is that the 
encoder is now complex: c = uG = u [ 1 1 P] (imagine if G 
was 9000 x 10000 so that P was 9000 x 1000) 

• more recently, Shu Lin and others have been looking at 
algebraic and combinatoric designs 

• we have been examining both approaches for CCSDS near- 
earth standards work 
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E K /N n in dB 

D 0 


A rate r = 0.82 MacKay LDPC code with BPSK7QPSK 
signaling on an AWGN channel, with and without RFI bursts 






Probability of Error 
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A rate r = 7/8 (= 0.875) quasi-cyclic LDPC code based on 
difference sets for BPSK/QPSK signaling on an AWGN 
channel, with and without RFI bursts. 
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Future Research 
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Therefore, the communications channel has 
limited bandwidth and may be error prone 


Introduction 
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The motion compensation process makes the 
complexity of their encoders considerably 
higher than that of their decoders 
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In-Subband Motion Compensation 
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Out-of-Loop Motion Compensation 



(a) is encoder; (b) is decoder 
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transform (e.g., wavelet) coefficient encoding 

• Heterogeneous communications networks 

• Channels with time-varying capacity 



Out-of-Loop Motion Compensation 
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Out-of-Loop Motion Compensation 
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Frame-Rate Scalable Motion 
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(a) is encoder; (b) is decoder 
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Bit rate: 9.2 % lower 

Mean Squared Error (MSE) in reconstruction: 29.6% 
lower 


Results 





0.39 bpp, MSE= 382 
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Developing new invertible motion 
compensation operators 
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Receivers: Why Implement Digitally? 
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Receivers: Why Parallelize? 
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Serial Implementation of QAM Receiver 
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Parallel Implementation of QAM Receiver 
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Possibility #1: Parallel QAM Receiver 
Based on Filterbanks (JPL, 1995) 
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Possibility #1: Parallel QAM Receiver 
Based on Filterbanks (cont.) 
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Possibility #2: Parallel QAM Receiver 
Based on Frequency Domain Filtering 

(GSFC, 1998) 
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IC Implementation 
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Review of Adaptive DSSS Receivers 
(Fractionally Spaced) 
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Review of Adaptive DSSS Receivers 
(cont.) 
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Parallel Implementation of Adaptive 
DSSS Receiver: 1st Approach (trivial) 
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Parallel Implementation of Adaptive 
DSSS Receiver: 2nd Approach 
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Parallel Implementation of Adaptive 











Parallel Implementation of Adaptive 
DSSS Receiver, 2nd Approach (cont.) 
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Parallel Implementation of Adaptive 
DSSS Receiver, 2nd Approach (cont.) 













Comparison of Serial and Parallel 
Receiver by Simulation using NLMS 
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Comparison of Serial and Parallel Receiver 
by Simulation using NLMS (cont.) 
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emarks to NLMS Simulation 
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Gain factors speed up filter convergence 

DCT, Hadamard transforms perform well, DFT poor 



Comparison of Serial and Parallel 
Receiver by Simulation using RLS 



BER MSE of Serial Receiver 





Comparison of Serial and Parallel 
Receiver by Simulation using RLS (cont.) 
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Comparison of Serial and Parallel 
Receiver by Simulation using RLS (cont.) 
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Remarks to RLS Simulation 
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Optimal Transform Design: General 
Approach 
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Optimization of Transform, T for 
Maximum Output SNR (cont.) 
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Performance Comparison of Parallel 
Receiver Using Different Transforms 
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Optimization of Transform for Fast 
Convergence: Single User Case 


cn 


Ch t rl 
0,0^ 

3 c ^ 

.a .2 - 
§ ts + 

'H N ^ 

Ch ^ 


+ 

CO Tt 

Co 


rf / N 

^ ^ 2J 


/I'k H" 

D O 

O Jr <n 

S .« 

£? £ 2 

§5 £ £ 

> o + 

C/5 V-> . . 

<2 c s 


CN 

>5 


3 -I 

£ £ 


a T3 m p 

o Cd H 

•rt <D «Jh 

i&i 

■I SJ Hi 


is also in numerator of equation for the SNR 





vr 

Vi? 

VJ? 

VJ? 

Vi^ 

VT 

v*? 

w 

VI? 

VI? 

v^ 

Vp 

vP 

W 

1U? 

VP 

Vtf 

VP 

vi? 

vP 

VP 

vp 


VP 

VP 

vP 

VP 


C/D 

a) 

c/3 

<3 

u 

to 2 

c 3 ^ 


tin 


c 3 

N 


O 

C/3 

S 3 

C 

CO 

#N 

5-h 

<D 

CW 

S 3 



Oh 

o 


a) 

o 

c 

o 

bO 

«-H 

a> 

> 

a 

o 

U 



VI? 


m 





Berner, De Leon: Parallel Digital Architectures for DSSS Receivers 


Computational Costs 
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- Synchronization issues, non-square transforms, more 
realistic channel conditions 
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TWT Amplifier Model 






Input voltage (normalized) 





AM/PM Characteristics of a Commercial TWT 



Input voltage (normalized) 



Spectra of Conventional RRC Filtered QPSK 
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FQPSK Block Diagram 










Waveform Set of EFQPSK 





Waveform Set of EFQPSK contd. 



s 7 (t) = sin(7 rt/T) —T/2 < t < T/2, s 15 (t) 



Vector Plot of EFQPSK 



Inphase 







Spectra of EFQPSK 
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Spectra of EFQPSK with Unbalanced Data 
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Table 1: Line strengths with unbalanced data. The reference spec- 
trum level is zero frequency for balanced data. 

















Spectral lines are produced when unbalanced data are used. 





Constrained Envelope Root Nyquist (CERN) 
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On-time samples: samples at time instants t = iT. 
Off-time samples: samples at time instants t — (i + 0.5 )T. 



CERN Generation 



Threshold 
sample value 












Magnitude Probability Density Function (8 PSK) 



Sample magnitude 



Magnitude Probability Distribution Function 
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QPSK CERN Spectra with Ideal Clipper 
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QPSK CERN Spectra with a Commercial TWT 
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Effect of the Predistorter on AM/ AM 



Input voltage (normalized) 








QPSK CERN Spectra with Predistorter 



f/Rb 




Conventional RRC QPSK Spectra with Predistort 
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Effect of sampling rate on EFQPSK 



Sampling rate (i 



BER performance of uncoded EFQPSK 



SNR (dB) 








BER Analysis/Simulation for uncoded EFQPSK 



Wiener (simulation) 
Average MF (simulation) 



Sensitivity to predistorter phase error (EFQPSK) 
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BER Performance of Coded EFQPSK) 



SNR (dB) 





Convergence of Adaptive Receivers for EFQPSK ) 



number of iterations 






BER Performance of QPSK CERN 
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improvement over CERN. 







Conclusions 
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beveral tuture research directions in terms ot improving th( 
performance of EFQPSK and CERN techniques have been 
indicated. 
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TWT Characteristic 
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Proposed Solutions 
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- Faster convergence than direct adaptation 

- Run TWT in linear region during adaptation 
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LMS and RLS performance nearly identical. 







Channel Modeling 
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Keep it as simple as possible, 
not simpler. 
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Keview proposed compensation techniq 
Develop or borrow simulation software. 
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First step is good channel modeling and estimation 
procedures. 



Lightweight Optical Wavelength 
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consumption: < ouum 

Day / night operation 

LEO communication 
without a laser in space 
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LOWCAL: Lightweight Optical Wavelength 
Communication without A Laser in space 





Locked: 7 m 





Comparison With Othe 
I Experiments 
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DCPK COMMUNICATIONS FORMAT 
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where c represents the speed of light 


DOWNLINK LINK EQUATION 



where: 

■^■scintillation represents the margin required to compensate for beam scintillation. 
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D retro represents the retro-modulator diameter 
0 represents the spacecrafts angle with Zenith 


DCPK SNR 
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k l represents tne load resistor. 

F t represents the noise figure of the amplifier, 
e represents the extinction ratio of the LCS 
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DATA RATE VERSUS POWER 
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Maximum angle from Zenith 60 degrees 
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DOWNLINK SUMMARY 
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With faster modulators data rates of Mb/s are possible 
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LIGHTWIRE FLIGHT SUBSYSTEMS 
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filter 
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LIGHTWIRE FORMAT 
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where c represents the speed of light 




Conceptual Block Diagram for Uplink Modulation 
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Lightwire FSK Demodulation 
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FSK Demodulation Theoretical 
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Phase 
Comparator 
Output 
(3 states) 


FSK Demodulation Experimental 
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Ij represents the intensity incident upon the spacecraft, 
m represents the modulation index. 
c°fsk represents the subcarrier frequency. 


UPLINK SNR CALCULATIONS 
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k represents Boltzmans constant. 

T represents the temperature in degrees Kelvin. 
R l represents the load resistor. 

F t represents the noise figure of the amplifier. 


PRELIMINARY UPLINK 
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Summary of Unique Features 
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Patent pending 






Lightwire: FSK and Circular Polarization Keying 

transparent link pair 



Receiver p I FAD OF ^ r*) \/ Polarizing Beam Splitter 



Subsystem Tests 
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Preliminary System Tests 
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Lightwire Demonstration 

(Full Duplex on a single beam) 

• BER < 1 .67 10~ 7 (no errors in 10 minutes at 10 kbps) 




errors were detected over a 10 minutes period for any 
of the above link tests 



GOALS FOR FY 01-02 
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Safety review preparations 





• Funding 2.5 Gbps link at 1 .5 microns 

• Funding a second NASA type link 

DOE (LANL) 

- Considering funding an additional experiment 



